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INTRODUCTION 



A rigorous ray tracing theory has already been developed to first order 

where it was shown that an effective sound velocity profile was in general 

dependent upon frequency.^ A higher order theory for wave mechanics has 

2 

also been developed in one dimension. Exact albeit numerical wave functions 
have recently been developed for bound states (i.e,, trapped modes in under- 
water acoustics) in one dimension where the turning point (vertex depth in 
underwater acoustics) of the phase integral have been shown to recede to >» 
for indices of refraction that remain finite and where the effective action 
variable (the phase integral around the branch cut between the turning points 

for the N^*^ eigenvalue, N = 1,2,3, ...) in units of phase has been shown to 

t h 3 

be N cycles for the N^" trapped mode. The modes have been shown to be 
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numerically exact within the limits of computer round-off. ’ 

In this exposition, the terminology "rigorous ray tracing" denotes an 
exact, albeit numerical ray tracing while "classical ray tracing" denotes 
the conventional ray tracing that is valid in the short wavelength or 
infinite frequency limit. 

Herein we apply some of these developments in wave mechanics to under- 
water acoustics in order to generate for depth-dependent sound velocity pro- 
files a rigorous ray theory that is exact, albeit numerical, and that is 
mutually consistent with a normal mode theory. The objectives of this inves- 
tigation are to gain insight into the physics of rigorous ray tracing from 
a few examples of the exact effective sound velocity profile, to develop the 
correct equations of motion for the ray, and to show that the ray directions in 
general are not orthonormal to the wave surfaces of constant .phase (i.e.. 
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the wave normals and the ray directions for the same points are not in 
general paral Ie1 ) . 

In this exposition we assume that the sound velocity profile is a func- 
tion of depth only and that for finite depths, the inverse square of the 
sound velocity profile is finite. The surface is at zero depth. The earth 
is assumed to be flat and we assume a bottomless ocean. Herein cylindrical 
coordinates are used and we assume the source is located on the z-axis at 
some depth, z^. Hence we have azimuthal symmetry. We also assume that the 
sound source is a point harmonic source whose frequency is ca in radians per 
second. All numerical computations were based upon a Runge-Kutta fourth- 
order technique with a step size of 25 m. 

In Section I of this exposition, we develop Hamilton's principal func- 
tion, S, for rigorous ray tracing which is based upon the alternative de- 
velopment presented in Appendix A of Ref. 1. (Heretofore S had been cited 
as the effective Hamilton's principal function.) Subsequently we develop 
the equations of motion by the Hamilton-Jacobi process and we develop the 
relationship between the ray direction and wave normal. In Section II we 
generate some sample exact effective sound velocity profiles, C^Cz) from 
a simple sound velocity profile, C(z). From these exact samples, we glean 
some insight into the behavior and application of rigorous ray tracing. A 
brief outline of applicable portions of classical ray tracing is presented 
in Appendix A. 
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I . THEORY 



Here we shall develop the equations of motions for a ray and show that 
the ray is in general not normal to the acoustical wave surfaces of constant 
phase. Here for convenience and without loss of generality, we express the 
sound velocity profile, C, which has been assumed to be a function of depth 
only, by the form 

C(z) = [A - 

where z is depth and where A, a real constant, and B(z), a real function, 
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are chosen such that at the source depth, z^, of the sound, C(z^) = A . 

As C(z) is only depth dependent, the Helmholtz equation is separable 
such that the z component may be expressed 

3 ^e(z)/ 8 z^ + J[k - C'^ - B(z)]0(z) = 0 (1) 

_2 

where e is a scalar, uj is frequency in radians per second, is the real 

_2 

separation constant such that A > C . In addition, C is also the vertex 
(WKB turning point) velocity. Let us use an ansatz, which has been inspired 
by the WKB approximation, for 6 of the form 

®± " (A - C'^ - exp[xi/ o3(A - C'^ - (2) 



Substituting Eq. (2) into Eq. (1) generates a nonlinear differential equation 
for our new variable, 3, expressed by 



_L 8^B/8z . 

4u)^ (A-C-2- 3) 



5 / 93/3Z 

IBoj^'A- C'^ - 3^ 
m 



B(z). 



(3) 
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The initial conditions for a require that Huygen's construction may be 
accomplished in the neighborhood near the source J Therefore the initial 
conditions are 



and 



e| = B(z ) = 0 

z=z 

s 



33/32 1 = dB/dzj 

z=z z=z 

s s 

From the form of Eq. (3), we see that 3 is a function of z, o and C^. We 

note that for trapped modes Eq. (3) has a nodal-type critical point^ such 
_2 

that 3 A - C from below as z ->- <». This node at z = <» moves the vertex 
m 

depth (turning point) for the ray (for flat earth) for rigorous ray tracing 
from the WKB vertex depth to infinity. Further investigation of 3 is de- 
ferred to Section II , 

The phase of our ansatz is given by 

/ uj[A - 3(2',u),C^) - 

= / a)[C’^(z'.o).C^) - C’^]^/^dz* = ooW(z,u).C^) (4) 

where by analogy to Eq. (A4) W is Hamilton's characteristic function for 
rigorous ray tracing and is defined to be the effective velocity given 
by 



Ce(z.u).Cm) 



[A - 3(z,cj,C^)] 



- 1/2 



(5) 
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We develop flamilton's principal function, S, for rigorous ray tracing as 



S = W(2,U),C ) + ^ (6) 

m 

because range is a cyclic coordinate and azimuth is symmetric. By Eqs. (4) 
and (6) and by using 8S/3r = (VC^) to change temporarily from a function 
of to 8S/8r, we have the Hami 1 ton- Jacobi equation for rigorous ray 
tracing expressed by 

9S/9r - (C'^(z,u),9S/9r) - (9S/9z)^]^^^ = 0. (7) 

There is a subtle difference between Eqs. (7) and (A2) because the range 

dependence, which only appears explicitly in Eq. (A2), appears in Eq. (7) 

_2 

both explicitly and in implicitly. Thus, even though is a separation 
constant for the Helmholtz equation, is no longer a separation constant 
for Eq. (7) (albeit Eq. (6) is still valid for cyclic r) , but still 

is a transformed momentum of the Hami 1 ton- Jacobi process for both Eqs. (7) 
and (A3). Hence is in general anisotropic. 

Either S or W (at the worker's option) may act as the generating func- 
tion from which the equations of motion for rigorous ray tracing may be 
derived. But we note that the ansatz may be expressed as 

9+ = (9W/9z)"^'^^ exp ±io)W. 

As the ansatz describes the normal mode (the analogous exact, albeit numeri- 
cal, quantum mechanical eigenfunctions have been so described elsewhere ’ ), 
we may derive the eigenfunctions from W by quantizing the product of to and 
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the action variable for rigorous ray tracing to be N cycles of phase, 

N = 1,2,3, ... . Hence W is the keystone by which we may derive mutually 
consistent normal mode and rigorous ray tracing theories thus rendering 
a unifying principle. 

The wave normals for rigorous ray tracing are given by the gradients 
of S which for cylindrical symmetry renders 



VS = 






c 1 . 

m r 



Thus the wave normal is embedded in a radial -vertical plane and the wave 
normal may be described by the deflection/elevation (D/E) angle, i}/. This 
D/E angle is given by 

ij; = tan'^(^)^- 1]^/^ = cos'U-^) (8) 

^e m 

which may be expressed in a manner similar to Snell's law for classical 
ray tracing as 



C 



m 



cos l|^(z) 



( 8 ') 



where is an implicit function as opposed to the classical ray tracing, 
Eq. (Al), where it is an explicit function. 

We now derive the continuous motion for rigorous ray tracing for the 
ray passing depth z. By the Hamil ton- Jacobi process the range r that a 
ray passes depth z is given for cylindrical symmetry by 



as 



0 



aw 

KT7y 



r 
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or 



r 



z 



s 




<lz 



(9) 



where the Hamil ton- Jacobi constant coordinate, which normalizes range with 

respect to source location, is zero for sources located on the z-axis. It 

is noted that the integral in Eq. (9) is a line integral and that one must 

integrate along the branch cut and around any encountered branch points 

(i.e., vertex depths) onto the next Riemann sheet. The node in Eq. (3) 

at z = oo for trapped rays induces C -> C from below and dC/dC -> 1 from 

e m e m 

below as z 00. Thus the vertex depth for the ray is at infinite depth. 
Nevertheless, the range quadrature, Eq. (9), for rigorous ray tracing 
remains bound and integrable for trapped rays as shown in Section II. 
Furthermore in the limit by L'Hospital 's rule, the integrand in 

Eq. (9) goes to zero. 

For cylindrical symmetry by Eq. (9), the ray also is embedded in a 
radial-vertical plane. The D/E angle, 4>,_for the ray for rigorous ray 
tracing by Eq. (9) is given by 



which in general differs from as given by Eq. (8). Hence the ray direc- 
tion is not parallel to the wave normal. There is a precedent for rays not 
being normal to the wave surfaces of constant phase or W in the optics of 






CO t { 



C,n-(vc^)’(av8gi^ 

2 2 1/2 ^ 



m e 



( 10 ) 
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anisotropic crystals. In a proceeding paragraph of this section, it has 
been shown that in general is anisotropic. 

At the source, 80^/30^^ = 0 in order to satisfy the boundary condi- 
tions for Eq. (3). Hence at the source, we have by Eqs. (8) and (10) 
that. = 4 >(z^) . 

Equation (10) is Snell's law for the ray for rigorous ray tracing. 

Because Eq. (3) has a node at z = we have that ->• 0 and 4> * 90° as 
z ->■ oo. Hence the ray is parallel to the wave surface at the vertex point 
at z = <» for rigorous ray tracing. 

It is noted that Ref. 1 needs correction because Eq. (8) of this exposi- 
tion was given in Ref. 1 as the applicable modification of Snell's law without 
distinguishing between wave normals and ray directions. Nevertheless all 
applications in Ref. 1 of Snell's law, which were based upon an effective 
sound velocity profile developed only to first order, are correct because 
to first order only, C is a function of both z and d but not C . Thus a 
first-order is isotropic. We note that if is independent of 
then Eq. (10) reduces to Eq. (8). 

We offer a caution with regard to rigorous ray tracing with an effec- 

- 1/2 

tive sound velocity profile as given by Cg«f*(A-8) where 8 is the solution 
to Eq. (3). The worker is admonished to use only in a Hamilton- Jacobi 
process (regrettably, it does not render correct results when is substi- 
tuted for C into a Lagrangian or Hamiltonian calculus of variation process). 
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II. ANALYSES OF SAHPLE APPLICATIONS 



Let us assume a sound velocity profile, C(z), such that 

C(z) = [(1420)^ - iL . Q 36 z j-1/2^ ^ ^ ^ 

{lAZO)"^ 

where C is in meters per second and z is in meters. The source depth is 

assumed to be zero. This profile was chosen for investigation because 

(1) it represents the profile that is dominated by the pressure gradient 

in the absence of a thermocline and (2) in a perturbation expansion the 

first order calculation of C is identical to C. Here we examine C for 

e e 

various values of and cj. Subsequently we rigorously trace a particular 
ray to exhibit its penetration of the classically forbidden region beyond 
the WKB vertex depth and to exhibit the divergence between the ray direc- 
tion and wave normal . 

Let us now numerically solve Eq. (3) and exhibit by Fig. 1 the effec- 
tive sound velocity profile, C^, for a harmonic surface source (z^ = 0) 
with 0 ) = 50 rad/s for rays with initial D/E angles [note il;(z^) = 4>(z^); 
z^ = 0] of 5®, 10®, 15® (i.e., = 1425.4242 mps , 1441 .9058 mps, 1470.0922 

mps respectively). In Fig. 1, the nodal behavior of Eq. (3) is manifested 
by -*■ C from below as z -► «>. Hence we observe that the vertex depth 
recedes to z = Well above the WKB vertex depths, and C are good 
approximations of each other (this may not be true for less-well-behaved 
profiles^), and by the initial conditions is tangent to C at the source 
depth z^. For completeness, the actual sound velocity profile, C, is also 
exhibited on Fig. 1. Note that in the limit where the initial B/E angle 
goes to 90®, (i.e., -> °°) , then would mimic C as closely as possible 

for the value of A (here A^^^ = 1420 m/s) but still would not go to C. 
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However in the limit A >• and the initial D/E angle is not 0°, then ^ C 
(the analogy in quantum mechanics would be the infinite energy case). 

In Fig. 2 we exhibit various effective sound velocity profiles for 
surface sources with an initial D/E angle of 10° (i.e., = 1441.9058 mps) 

and with oj = 5 rad/s, 50 rad/s, 500 rad/s. The sound velocity profile, C, 
is also exhibited on Fig. 2. We note that for finite depths above the 
WKB vertex depth approaches C as u) increases and that for finite depths 
below the WKB vertex depth, approaches as oj increases. In the limit 
cj 00 , we have that 



Lim C 

e 

UJ-KO 




c < C„ 
— m 




C > C 



m 



( 11 ) 



Let us now investigate the correspondence principle of how rigorous 
ray tracing goes to classical ray tracing in the high frequency limit. 

In the case oj -► «> and below the WKB vertex depth (which is the classical 
ray tracing vertex depth) in the classical forbidden region, the ray direc- 
tion is parallel to the wave front and so the ray is embedded in a wave 
surface of constant phase. Hence below the WKB vertex point for oo ^ 
the ray senses a constant phase and thus the ansatz, Eq. (2), in the forbidden 
region has a constant phase despite oj ^ that is, the integrand of Eq. 

(4) behaves such that (A - C~ - 3) -► 0 as oj The boundary condition 

for infinite frequency for the solution to Eq. (1) requires that in the 
classical forbidden region that e(z) = 0 for z below the WKB vertex depth. 
Therefore we have that 



9 = (0^ - 0_)/2, OJ -> «> 
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to I a i..‘low t.iu‘ .,1..*: Vv*:t x 



which goes to zerc^by L'Hospital's rule despite the singularity in the 

amplitude factor in both 0^ and 9 as given by Eq. (2). Thus for u -*■ 

there is no penetration below the WKB vertex depth, and by Eq. (11) below 

the WKB vertex depth the integrand of the range quadrature, Eq . (9), for 

rigorous ray tracing is zero by L'Hospital's rule while above the WKB 

vertex point Eq. (9) reduces to the range quadrature, Eq. (A5), for classical ray 

tracing as 9^^790^^ -> 0. Thus in the limit to -> «>, rigorous ray tracing corres- 

ponds to classical ray tracing [the quantum mechanical analogy is the corres- 

2 

pondence principle for (h /2y) 0 where fS is Planck's constant and p is mass]. 

We exhibit on Fig. 3 the rigorous ray trace for a surface source for 
to = 50 rad/s and the cj)(z^) = = 10° (i.e., = 1441.9058 mps) by 

numerically evaluating Eq. (9). The corresponding classical ray trace is 
exhibited also on Fig. 3 for a half-cycle in range (i.e., until the first 
WKB vertex depth is encountered at z = 1189.1117 m) where classically 
r = 13,490.938 m. While the corresponding rigorous ray trace has its 
vertex depth at z = <» which is not exhibited, we report that the range for 
the corresponding rigorous ray trace to the first vertex depth has been 
evaluated to be 13,557 ± 6 m. This infers that classical ray tracing for 
one-half cycle for these conditions has an estimated error of 65 m in range 
or about 0.498% in range. 

We note that while the vertex depth for rigorous ray tracing has re- 
ceded to infinity, the range quadrature has remained finite. We also report 
that the time quadrature, which is the action variable for a full cycle 
remains finite as inferred by the analogous quantum mechanical effective 

action variable remaining finite for bound states of the symmetric bi-linear 

3 4 

potential whose eigenfunctions are offset Airy functions. ’ 

On Fig. 4 we exhibit as a function of depth the wave noniial, \p, and 
the ray angle, <j>, for the ray that was rigorously traced on Fig. 3. For 
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classical ray tracing, which also is exhibited on Fig. 4, the wave normal 
and ray angle are equivalent. We note that well above the WKB vertex depth 
in the classically allowed region both and 4) approach their classical ray- 
tracing prediction. But below the WKB vertex depth we have that as z 
then '4; -> 0 and (}) -> 90°. Hence, as z the ray travels parallel to the 

wave front as noted earlier. 
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APPENDIX A 



CLASSICAL RAY TRACING 

In this appendix, we outline the development of the Hami 1 ton- Jacobi 
theory of classical ray tracing in the short wavelength. A, or high frequency 
limit (i.e., A ^ 0 or cd -»■ <») and show that the ray path is always normal to 
the wave fronts. Here we assume that the sound velocity profile, C, is 
only dependent upon depth, z, that the earth is flat and that the ocean 
is bottomless. We use cylindrical coordinates and assume the source is 
located on the z-axis at depth z^ . 

We start with the Lagrangian, L, for ray tracing that for azimuthal 
symmetry is given in cyclindrical coordinates by 



L(z,i) 



0 . 

crn 



where z = dz/dr where r is the radial coordinate and for sources on the z- 
axis r determines range. The conjugate momentum, P^, with respect to . 
depth is given by 



P 



z 




C(z)(l + 



sin 4) 

“cM 



where 4> is the depression-elevation (D/E) angle of the ray (for azimuthal 
symmetry, the ray is embedded in the radial -vertical plane). As L is not 
an explicit function of range, r, the classical ray-tracing Hamiltonian, H, 
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is given by 






cos (}) _ 1 

-cuf ■ -c^ 



(Al) 



which gives Snell's law since H is a constant of the motion for H indepen- 
dent of r. The constant of the motion for any particular ray may be 
expressed in terms of its vertex velocity, C^. The Hami 1 ton- Jacobi equa- 
tion for classical ray tracing is given by 



aS/3r + H(z,3S/3z) = 0 



or 



3S/3r - [C(z)‘^ - (3S/3z)^]^/^ = 0 (A2) 

where S is Hamilton's principal function for classical ray tracing and has 
the dimensionality of time. Thus the surface of constant S represents a 
wave surface. Because the only range dependence is restricted to the 
3S/3r term in Eq. (A2), we may assume a solution of the form 

S(z,C^,r) = W(z,C^) + r/C^ (A3) 

where W is Hamilton's characteristic function for classical ray tracing. 
From Eqs . (A2) and (A3), W may be expressed as 



14 



(A4) 



w M [c 



-2 






dz , 



The wave normal is given by the direction of the gradient of Hamilton's 
principal function, S, which by Eqs. (Al) and (A2) is given for azimuthal 
symmetry as 



VS 




r_L . _L,V2 
I 2 2 J 

m 



For azimuthal symmetry the wave normal is embedded in a radial -vertical 
plane and the wave-normal may be described by the D/E angle, of the 
wave normal as given by 



tan 








But the equation of motion for the ray for classical ray tracing from a 
source on the z-axis is given for azimuthal symmetry by 



3S 

alTTy 



0 




-C dz' 
[C^ - 



+ r 



(A5) 



where the z integral is a line integral with branch points. For azimuthal 
symmetry, the ray is also embedded in the radial -vertical plane and the 
ray's D/E angle, (p, is given by 

(p = tan'\dr/dz)”^ = tan[-y- 1]^^^ = 

r 
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Hence for classical ray tracing, the 
wave normal are equal to each other, 
the ray path is always normal to the 



D/E angle for both the ray and the 
Therefore by azimuthal synmetry, 
wave fronts. 
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FIGURE CAPTIONS 



1. The effective sound velocity profiles for a harmonic source at zero 

depth with w = 50 rad/s for initial D/E angles of 5°, 10° and 15° 
(i.e., = 1425.4242 mps, 1441.9058 mps and 1470.0922 mps respec- 

tively). The sound velocity profile, C, is represented by the 
dashed curve. 

2. The effective sound velocity profiles for a harmonic source at zero 

depth with an initial D/E angle of 10° (i.e., = 1441.9058 mps) 

for u) = 5 rad/s, 50 rad/s and 500 rad/s. The sound velocity profile, 
C, is represented by the dashed curve. 

3. The rigorous ray trace for a harmonic sound source at zero depth, with 
( 1 ) = 50 rad/sec, with an initial D/E angle of 10° for the sound 
velocity, C, exhibited on Figs. 1 and 2. The classical ray trace is 
represented by the dashed curve. 

4. The deflection/elevation angle for the ray, (j), and wave normal, ip, 
as a function of depth for the rigorous ray trace exhibited by 
Fig. 3. The deflection/elevation angle for classical ray trace is 
exhibited by the dashed curve. 
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